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A new approach to the analysis of three body decays is presented. Model-independent results 
are obtained for the S-wave Kir amplitude as a function of Kn invariant mass. These are 
compared with results from K~n^ elastic scattering, and the prediction of the Watson theo- 
rem, that the phase behavour be the same below Kr]' threshold, is tested. Contributions from 
1=7; and I — ^ are not resolved in this study. If 7 = ^ dominates, however, the Watson 



2 t^iiv_i ^ — I are not resolved in this study, 
theorem does not describe these data well. 



1 Introduction 



Decays of heavy-quark mesons are regarded as a potential source of information on tlie liglit- 
quark mesons they produce. For decays to three pseudo-scalar final states, kinematics and 
angular momentum conservation favor production of S-wave systems, so improving our knowl- 
edge of the particularly confusing scalar meson (J^ = 0"*") spectrum may be possible when 
the anticipated large, clean samples of such decays of D mesons from the B factories and the 
Tevatron collider become available. Extracting this information has, however, been done in 
model-dependent ways that make assumptions about the scalar states observed. Such assump- 
tions can influence the results, so new approaches are required. 

In this paper, we present a model-independent approach to the S-wave system in a study 
of the decays K~t:~^tt^ " observed in data from Fermilab experiment E791. We compare 

the S-wave amplitudes so obtained with our earlier, isobar model analysi^and also with K~-k~^ 
scattering. 

Measurements of K~it^ scattering come principally from SLAC experiment E135 (LASS)^, 
and cover the invariant mass range only above 825 MeV/c^. Data exist below this range, but 

"Throughout the paper, charge conjugate states are implied unless explicitly stated otherwise. 



with less precisionl^. More information in the low mass region is required if the possibility of 
the existence of a k state is to be properly evaluated. 

2 Data Sample 

The selection process for events used in this analysis is described in Refffl A signal consisting of 
15,079 — > i^~7r+7r^ decays, with a purity of ~ 94%, is obtained. Fig. ^shows the Dalitz plot 
with K~t:~^ squared invariant mass s plotted vs. s' . Horizontal (and the symmetrized vertical) 
bands corresponding to the i<'*(892) resonance are clearly seen. A striking and complex pattern 
of both constructive and destructive interference is seen near 2 (GeV/c^)^ due to either i^'Q(1430), 
i^J(1410) or Er|(1430). There is also evidence for Erj*(1680), difficult to see due to smearing of 
the Dalitz plot boundary resulting from the finite resolution in the three-body mass. 




Figure 1: Dalitz plot for D"*" — > i('~7r^7r^ decays. The squared invariant mass s is plotted against s' . The plot 
is symmetrized, each event appearing twice. Lines in both directions indicate values equally spaced in squared 
effective mass at each of which the S'-wave amplitude is determined by the model-independent partial wave analysis 
(MIPWA) described in section |3 Kinematic boundaries for the Dalitz plot are drawn for three-body mass values 
M = 1.810 and M = 1.890 GeV/c^, between which data are selected for the fits. 

The most striking effect observed is the asymmetry in the if* (892) bands, most easily 
described by interference with a significant S'-wave contribution to the decay. In this paper 
we are able to extract information on the S-wave using the -fC*(892), and also the other well 
established resonances in the Dalitz plot, as an interferometer. 



3 Method 

In Ref.^, as in most earlier analyses of D decays to three pseudo-scalar particles ijk, we use the 
"isobar model" . Details of this are given in^^. In this model, the decay amplitude A is described 
by a sum of quasi two-body terms D — > R+k, R — > i+j, in each of the three channels k = 1,2, 3: 



^ FR{p,rR, J) 
mt — Sii — irrir 



A = doe''' + E due''" — ; , , X F,(g, r„ J) Mj{p, q) (1) 

n=l 



In this, Sij is the squared invariant mass of the ij system. J is the spin, the mass and 
^Rni^ij) the width of each of the resonances Rn seen to be contributing to the decay. Fr 



and are form factors, with effective radius parameters r^j and r^,, for all Rn and for the 
parent D meson, respectively, p and q are momenta of i and k, respectively, in the ij rest 
frame. Mj{p,q) is a factor introduced to describe spin conservation in the decay. The complex 
coefficients dnC^^" (n = 0, A^) are determined by the D decay dynamics and are parameters 
estimated by a fit to the data. The first, non-resonant (NR) term describes direct decay to 
i + j + k with no intermediate resonance, and do and are assumed to be independent of Sij. 
For — > K~tt^7t'^ decays we Bose-symmetrize A with respect to interchange of vr^ and tt^. 

In Ref.^we reported that the NR term was smaller than previously thought, and that a fur- 
ther term, parametrized as a new J = resonance k(800) with = (797±19±43) MeV/c^and 

= (410 ± 43 ± 87) MeV/cP, gave a much better description of the data. Here, we examine the 
K~Tr^ S-wave in a model-independent way. The 5- wave part of Eq. ^ (all terms with J = 0, 
including the NR term) is factored 

S = SisKn) >^ M^{p,q)F^iq,r^) =lnterp {cke'^'')xM^{p,q)F^{q,ro) (2) 

into a partial wave S{sktt), describing K~7r^ scattering, and the product MQ{p,q)F,j{q^rj^) 
describing the D decay. S(sA-jr) is interpolated between a set of points c^e*'^''' defined at 40 
K~TT^ invariant mass squared values indicated by the lines in Fig. ^ Each Ck and 7fc is 
regarded as an independent parameter determined by the data. 
We factor the P- and D- reference waves in the same way: 

V = F{sK^)xM^{p,q)Fo{q,ro);V = -D{sKn)xM^{p,q)Fj,{q,rn), (3) 

however, we parametrize the partial waves P{sKn) and D{sKn) exactly as in Eq. Q 

We make an unbinned likelihood fit to the data, Using the method described in Ref.^. This 
incorporates an incoherent background function describing the 6% of our sample not correspond- 
ing to true D decays. We measure 86 parameters - all (0^,7^) and the coefficients d^e^^^ for 
is:*(892), Kl{lQm) in the P-wave and K|(1430) in the L»-wave. For the K*(892), we define 
df^e^^*" = 1 to provide the reference phase. 

The fit results in an excellent description of the data. Comparison of the observed and 
predicted population of the Dalitz plot gives a probability of 50% for 363 bins. 

4 Results 

The 5"-, P- and D-waves resulting from the fit are shown in Fig. [2 They are compared with the 
model-dependent fit from Ref. The main S'-wave features of both fits agree well. Resonant 
fractions and the total iS-wave fraction (about 75%) also agree within statistical limits. 

We turn now to a comparison of the S'-wave amplitudes S(sx7r) measured here with the 
amplitudes T{sktt) measured in A'~7r"^ elastic scattering. We expect, for each partial wave J 
(for each iso-spin /) that S(sA'7r) = \/sk-k/p^'^'^'^'^Q{sk-k)'~^{sk-k) where Q{sktt) describes the 
dependence of K~n^ production in D decays on sk-k- The Watson theorem!^ requires that, 
provided there is no re-scattering of the K~t:'^ from vr^, that Q is a real function, so that 
phases found in D decay should match those in K~tt~^ elastic scattering data. 

1=1/2 phases measured by LASS are plotted in Fig.|21 There is a large offset in the 5-wave, 
about 75°, not seen in P- or L)- waves. The shapes of S- and P- waves are also not the same. 
Unless significant admixture of / = 3/2 K~tt~^ production occurs, these results suggest that the 
conditions for the Watson theorem are not met in these data. 
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Figure 2: (a) Phases (left) and magnitudes (right) of S-wave amplitudes for the K~n'^ systems from _D+ 

tt+tt"*" decays obtained from the MIPWA fit described in the text. The effect of adding systematic uncertainties 
in quadrature with statistical errors is indicated by extensions on the error bars. The P-wave and D-wave 
amplitudes are plotted in (b) and (c), respectively, as curves derived from the isobar parameters used in these 
waves, and from the full error matrix resulting from the fit. Shaded areas represent one standard deviation limits 
on the amplitudes. In all plots, the dashed curves show one standard deviation limits for the prediction of the 
isobar model fit. 1 = 1/2 phase measurements from the LASS experiment, with vertical and horizontal error bars, 
are included in the phase plots. The vertical lines mark the K~rj' threshold, the upper limit of elastic scattering. 
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